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Abstract 

Many applications of simulation models and related decision support tools for agriculture and 
natural resource management require daily meteorological data as inputs. Availability and quality of such 
data, however, often constrain research and decision support activities that require use of these tools. 
Daily solar radiation (SRAD) data are especially problematic because the instruments require electronic 
integrators, accurate sensors are expensive, and calibration standards are seldom available. The Prediction 
Of Worldwide Energy Resources (NASA/POWER; power.larc.nasa.gov) project at the NASA Langley 
Research Center estimates daily solar radiation based on data that are derived from satellite observations 
of outgoing visible radiances and atmospheric parameters based upon satellite observations and 
assimilation models. The solar data are available for a global 1° x 1° coordinate grid. SRAD can also be 
estimated based on attenuation of extraterrestrial radiation (Q 0 ) using daily temperature and rainfall data 
to estimate the optical thickness of the atmosphere. This study compares daily solar radiation data from 


NASA/POWER (SRAD np ) with instrument readings from 295 stations (SRAD 0 b), as well as with values 
2 that were estimated with the WGENR solar generator. WGENR was used both with daily temperature and 
precipitation records from the stations reporting solar data and records from the NOAA Cooperative 
4 Observer Program (COOP), thus providing two additional sources of solar data, SRAD W g and SRAD C o- 
Values of SRAD NP for different grid cells consistently showed higher correlations (typically 0.85 to 0.95) 
6 with SRADob data than did SRAD W g or SRAD C o for sites within the corresponding cells. Mean values of 
SRADob, SRADwg and SRAD NP for sites within a grid cell usually were within 1 MJm" 2 d"' of each other, 
8 but NASA/POWER values averaged 1.1 MJm “d" 1 lower than SRADob- The magnitude of this bias was 
greater at lower latitudes and during summer months and may be at least partially explained by 
10 assumptions in ambient aerosol properties. Overall, the NASA/POWER solar radiation data are a 
promising resource for regional modeling studies where realistic accounting of historic variation is 
12 required. 

14 Abbreviations 

COOP, NOAA National Weather Service Cooperative Observer Program; NASA/POWER, NASA 
16 Prediction Of Worldwide Energy Resources; Q 0 , daily integral of extraterrestrial insolation; RMSE, root 
mean squared error; SRAD, daily integral of solar radiation; T max , daily maximum temperature; T min , daily 
1 8 minimum temperature 
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1. Introduction 


2 

Many agricultural and natural resource management efforts involve spatial scales above the field 
4 and fann levels. Applications range from monitoring regional water use, to identifying promising zones 
for production of new crops, to targeting of specific cultivars or crop traits, to determining the potential 
6 impact of climate change and potential options for adaptation. Spatial assessments often consider climatic 
variation and increasingly, long-term records of daily weather data are required to examine climatic risks 
8 or trends related to climate change. Such analyses, however, are usually constrained by the availability 
and quality of the observed long-tenn meteorological data. Weather stations may not be available in the 
10 regions of interest, and individual stations may lack data for long time intervals. Weather data per se may 
show local variation due to positioning of the station and the instrument, instrument calibration drift, 

12 change in instrumentation, and other factors (Younes et al., 2005; Davey and Pielke, 2005). Solar 
radiation data have long been recognized as especially problematic (Durrenberger and Brazel, 1976; 

14 Stoffel et al., 2000). Radiation must be correctly integrated at low sun elevation angles and over all 

wavelengths. Radiometers using thermopiles are expensive, while lower-cost silicon pyranometers are 
1 6 less accurate. Both types of sensors require electronic circuitry to integrate readings over time and are 
sensitive to ambient temperatures. Sensor calibration is difficult because accurate reference values 
1 8 (besides 0) cannot be produced through simple techniques; thus sensors are usually cross-calibrated to 
radiometers whose calibrations are traceable to standards such as those maintained by the National 
20 Institute of Standards and Technology. 

The Prediction Of Worldwide Energy Resources (NASA/POWER) project at the NASA Langley 
22 Research Center provides daily data for surface solar radiation and other weather variables on a 1° x 1° 
geographic coordinate grid for the entire globe (Stackhouse, 2010a; See Table 1 for an overview of the 
24 data available from the POWER archive.). The solar data are inferred from satellite observations of the 
outgoing top of atmosphere (TO A) radiances via an updated version of the Pinker and Laszlo (1992) 

26 radiative transfer based algorithm that was used to produce the fluxes for the Global Energy and Water 

Cycle Experiment (GEWEX) Surface Radiation Budget (SRB) project solar algorithm v2.81 (Gupta et al., 
28 2006). Within this algorithm, a calculated TOA albedo is matched to an inferred TOA albedo from 

measured dark clear-sky background and instantaneous (every 3 hours) clear-sky and cloudy-sky satellite 
30 visible radiances using a radiative transfer model (through the use of lookup tables) on a l°xl° degree 
grid. Using the background clear-sky radiance and information about the atmosphere (e.g., water vapor 
32 and ozone) the radiative transfer model infers an absolute surface albedo for a particular time and 

location. This step assumes a background aerosol and an assumed spectral albedo shape based upon the 
34 most prevalent surface type of the area. Then the surface albedo and the other input information are used 
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to infer the cloud and aerosol optical depths needed to match the observed clear and cloudy sky TOA 
albedos within a certain tolerance. The algorithm computes solar irradiance for clear and cloudy 
conditions using the inferred clear and cloud optical depths respectively and the other atmospheric 
information. The total solar irradiance is computed as the weighted sum of the clear and cloud fluxes and 
is integrated over the day to produce the daily averaged fluxes. The NASA/GEWEX SRB solar data 
currently in the POWER archive span the time period from July 1, 1983 through December 2007. 
Appended to this time series are irradiances estimate from the NASA CERES (Clouds and Earth Radiant 
Energy System) FLASHFlux (Fast Longwave and SFIortwave radiative Fluxes from CERES and 
MODIS) data sets (Stackhouse, 2010b) spanning from Jan 2008 through within one week of the current 
date (for algorithm description see, Kratz et al., 2010; the latter data are not considered in this paper). 
These data may be downloaded with other variables in a format compliant with the standards of the 
International Consortium for Agricultural Systems Applications (ICASA; www.ICASA.net , Hunt et al., 
200, 2006), which facilitates use in decision support tools such as the Decision Support System for 
Agrotechnology Transfer (DSSAT; Hoogenboom et al., 2010). 

Although initially developed for applications related to solar energy, energy consumption, and 
energy conservation, the NASA/POWER data appear suitable for agriculture and natural resource 
management (White et al., 2009; Bai et al., 2010). Their coarse geographic scale, however, may limit 
their usefulness: one degree of longitude is approximately 1 10 km at the equator and 80 km at 45° 
latitude. In assessing the effect of spatial resolution of precipitation and radiation data on regional yield 
forecasts, de Wit et al. (2005) concluded that a 50 x 50 km grid provided an adequate resolution. 
Similarly, for climate change research in the contiguous US, Janis et al. (2004) concluded that a network 
of 327 stations was adequate to monitor a 0.1 0°C decade' 1 temperature trend. Besides spatial resolution, 
there remains the question of whether the data assimilation process introduced important bias or other 
error in the data. 

Faced with a lack of reliable solar radiation data, numerous researchers have opted for generating 
values using data on latitude, temperature and precipitation as inputs. The procedures first estimate the 
daily extraterrestrial insolation (Q 0 ) based on latitude, date and the solar constant. This value is then 
reduced based on atmospheric transmittance or similar considerations (e.g., Richardson, 1981; Bristow 
and Campbell, 1984; Richardson and Wright, 1984; Hodges et al., 1985; Cooter and Dhakhwa, 1995; Liu 
and Scott, 2001). Transmittance is typically estimated from region-specific relations that consider the 
diurnal range of air temperature, which may further be varied depending on whether precipitation 
occurred that day. This paper compares the NASA/POWER solar radiation data with data from weather 
stations reporting instrument-based observations and from an implementation of the WGENR solar 
radiation generator (Garcia y Garcia and Hoogenboom, 2005; Garcia y Garcia et al., 2008). 
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2 2. Materials and Methods 

4 Daily data for solar radiation from NASA/POWER (SRAD NP ) were obtained from the web site 

(power.larc.nasa.gov; Stackhouse, 2010a), which allows for downloading of the data in several ASCII 
6 based formats. The dataset covered the continental US on a 1° x 1° latitude and longitude grid, 

representing 867 grid cells. The time interval considered was from 1 Jul. 1983 to 31 Dec. 2004. The 
8 NASA/GEWEX SRB solar version v2.81, as identified above, corresponds to the dataset used in our 
previous comparison of temperature data (White et al., 2008). 

10 Observed solar radiation data (SRAD 0 b) were obtained mainly from Internet sources such as state 

or regional climate networks (Supplement Table 1), which typically report data from automated weather 
12 stations using silicon pyranometers that output a current signal. Instantaneous values are registered and 
integrated digitally. For the widely used LI-200 Pyranometer 1 , LI-COR (2005) states that these sensors 
14 are calibrated against an Eppley Precision Spectral Pyranometer (PSP) using natural daylight, and the 
maximum absolute error is typically ± 3%, with a maximum of ± 5%. Datasets from stations were 
16 rejected if they provided less than two years of data between 1983 and 2004, in order to match the period 
represented in our NASA/POWER dataset. Data reported for several stations were clearly incorrect, 

18 including values much larger than Q 0 , negative values and values with a large systematic bias. Where 
detected, problem values were excluded based on the following criteria: SRAD 0 b greater than Q 0 , 

20 SRADob less than 0.2 MJm 2 d"', or time series of SRAD 0 b showing large, systematic deviations from 

patterns observed in other years, values of Q 0 or nearby locations. After the quality control process, a total 
22 of 295 stations were available, which were located in 181 grid cells of the NASA/POWER dataset (Fig. 

1 ). 

24 The WGENR program (Hodges et al., 1985; Garcia y Garcia and Hoogenboom, 2005), which 

uses the Richardson approach (Richardson and Wright, 1984) to estimate daily values of SRAD, was used 
26 to obtain two additional estimates of solar radiation from observed values for daily maximum and 

minimum temperatures and daily precipitation. The first estimate (SRAD W g) was obtained using daily 
28 temperature and precipitation records from the datasets of the observed values of solar radiation and thus 
coincided with the source locations of the SRADob- A parallel set of daily data (SRAD C o) were estimated 
30 from daily temperature and precipitation records from 855 individual ground stations from the National 
Weather Service Cooperative Observer Program (COOP). The COOP stations were selected based on 

1 Mention of a trademark, proprietary product, or vendor is for information only and does not constitute an 

endorsement by the USDA, NASA or the University of Georgia. 
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their being nearest to the centroid of a given grid cell of the NASA/POWER data set (White et ah, 2008). 
2 Use of the COOP observations provided SRAD C o for 855 grid cells of the NASA/POWER dataset. 

Comparisons of the solar radiation were based primarily on Pearson product-moment correlations 
4 calculated for daily data from individual stations, which were calculated using the Correlation procedure 
(PROC CORR) of the SAS 9.2 TS (SAS Institute Inc., Cary, NC, USA). Since annual variation in SRAD 
6 is dominated by readily calculated variation in Q 0 (e.g., Bristow and Campbell, 1984), the analyses also 
considered relations with Q 0 . 

8 

3. Results 

10 

3.1. Comparisons at a single location 

12 

A single, arbitrarily selected location, Immokalee, Florida (latitude 26.46°, elev. 1 1 m) was used 
14 to illustrate comparisons for one location. Based on correlations (Table 3), the SRAD NP showed the best 
agreement with SRAD 0B , with a correlation of 0.86 (P < 0.001). While the overall good agreement is 
16 borne out by Fig. 2, values of SRAD NP were consistently lower than for SRAD ()B , and the means of all 
daily values were 19.5 MJd ‘m 2 for SRAD 0B and 18.0 MJd 'rn 2 for SRAD NP . Mean values of SRAD WG 
18 and SRAD C o were 18.0 and 18.1 MJd'm' 2 , respectively. Values of SRAD 0B included 7 daily values out 
of 3151 that exceeded 95% of Q 0 , but excluding these values had minimal effect on the correlations and 
20 means. 

22 3.2. Comparisons over all sets of obser\>ed solar radiation data 

24 Across the 295 locations considered, SRAD NP exhibited higher correlation with daily variation in 

SRAD 0B , with many correlations of 0.9 or higher (Fig. 1 A and Table 3), than SRADwg or SRAD C o. 

26 Correlations between SRAD 0B and SRADwg were typically between 0.8 and 0.9 (Fig. IB and Table 3), 
while correlations between SRAD 0B and SRAD C o were slightly lower (Table 3). Interestingly, 

28 correlations of Qo with SRAD 0B were similar in value to those from the two weather generators (Table 3). 

Values of root mean square error (RMSE) for prediction of SRAD 0B by SRAD NP were generally 2 to 3 
30 MJm 2 d _1 (Fig. 3a) while RMSE values for SRAD 0B and SRAD C o were 4 to 5 MJm" 2 d"' (Fig. 3b and 3c). 

These results suggested that SRAD NP data represented day-to-day variation in SRAD 0B better than values 
32 from WGENR, a conclusion also supported by density plots comparing SRAD 0B with SRAD NP and 
SRAD C0 (Fig. 4). Thus, the data assimilation process used with the NASA/POWER data appeared 
34 superior to approaches that try to recreate variability in SRAD by considering daily maximum and 
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minimum temperatures and precipitation patterns (e.g., wet or dry days) as done in weather generators. As 
2 an aside, we note that Figure 4B also evidenced a problem with WGENR in that it appeared to produce an 
excess of values around 6 MJm 2 d _1 and which coincided with days with rainfall (data not shown). 

4 Comparisons of means indicated that the NASA/POWER values were slightly lower than 

observed values, with a mean across all stations of 16.2 for SRAD NP vs. 17.4 for SRAD 0 b (Table 2 and 
6 Fig. 1C). This difference was greatest in the summer months and was more pronounced at lower latitudes 
(Fig. 5A and 5B). However, when expressed on a relative basis (Fig. 5C), the differences were more 
8 pronounced in winter months. We note that SRAD WG and SRAD C o also tended to have mean values less 
that SRADob (Table 3 and Fig ID). 

10 A possible source of discrepancies in SRAD NP values relative to SRADob might relate to 

elevations of individual weather stations as compared to the elevation for NASA/POWER grid cell, which 
12 is the average of the topography associated with the 1 -degree cell. Since the thickness of the atmosphere 
decreases with elevation, clear-sky transmittance increases with elevation. The elevation of individual 
14 weather stations differed from the average elevation of the associated grid cell from the NASA/POWER 
dataset by as much as 800 m. Comparisons of correlations of SRAD 0 b with SRAD NP showed a weak 
16 trend related to differences in elevation (Fig. 6A). There was a slight relation between elevation 

difference and difference in mean values of SRAD 0 b and SRAD NP (Fig. 6B), but given that the mean 
18 elevation difference was only 219 m (Table 3), the net bias due to elevation differences would be less 0.4 
MJm‘ 2 d _1 . 

20 

3.3 Comparisons of NASA/POWER data with data generated using NOAA COOP data 

22 

The availability of the large set of data for paired NASA/POWER grid cells and NOAA COOP 
24 locations allowed a more detailed comparison for the continental US. The correlations between SRAD NP 
and SRADco were largest in the western USA and lowest in southeastern regions (Fig. 7A). The overall 
26 mean value of SRAD NP was 15.0 Mld^m 2 vs. 15.9 MJd'm 2 for SRAD C o (Table 2), again suggesting that 
values of SRAD NP are lower than other estimates. Mean values of individual cells diverged by as much as 
28 3 MJ d’Nrf 2 (Fig. 7B). The best agreement for means occurred in California, Oregon and southern to 

eastern states. In the Rocky Mountain region, mean values of SRADco were large relative to SRAD NP , 

30 while SRADco values were low for certain coastal regions. 

32 3.4 Data quality issues 
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In the initial analyses, several locations showed correlations between SRAD 0 b and SRAD NP that 
2 were less than 0.8, suggesting possible problems with the observed values. Data from Semmes, AL 
(Supplement Figure 1A) are indicative. For this location, SRAD 0 b showed a steady decline from 1995 
4 through 2004, so the location was excluded from the analyses. Plots comparing SRAD 0 b to Qo proved 
difficult to evaluate visually, so data were re-plotted as the ratio of SRAD 0 b to Qo (Supplement Figure 2). 
6 Although not analyzed quantitatively in this study, it appeared that the maximum value of this ratio, 

which would correspond to very dry, clear-sky conditions, is approximately 0.8 (as assumed in WGENR), 
8 and this value is marked by a reference line on the graphs. 

For seven weather stations, low correlations of SRAD 0 b with SRAD NP were not associated with 
1 0 errors visible in time series plots of SRAD OB , and the data showed mean values and patterns of variation 
similar to neighboring sites. Discrepancies in observation time are problematic in reporting of daily 
12 temperature data, so we tested whether correlations of SRAD 0 b with SRAD NP improved if it was assumed 
that the reported date of observation was a day later than the actual date. This assumption improved the 
14 correlations for these stations from a mean value 0.62 to a mean of 0.92 (Table 4), implying that the 

weather stations were reporting solar data with a one day offset (delay). 

16 

4. Discussion 

18 

The comparisons of the different sources of solar radiation data suggested that data from 
20 NASA/POWER reproduced variability at a daily time scale better than either set of generated values 
(Table 3 and Fig. 1). Comparisons of mean values (Tables 2 and 3 and Fig. 1), however, indicated that 
22 SRAD np values were often lower than observed values or values derived from WGENR. Figures 1, 5B, 

5C and 7 suggest that there is a regional component to the bias, although the variation in Fig. 7 also may 
24 result from bias in SRAD C o values. The subsequent version of the NASA/POWER solar irradiance data 
set showed reductions in this bias that are directly attributable to a reduction in the background aerosol 
26 specification. The new aerosol climatology based upon an upgraded NC AR Model of Atmospheric 
Transport and Chemistry (MATCH, Rasch et al., 1997) that resulted in the initial background aerosol 
28 optical depths being reduced from about 10% in the southeastern US to over 50% in the northwestern US. 

Although the background aerosol does not determine the final aerosol optical depth for a given space and 
30 time (as noted above), it governs the optimization of the match between the inferred and computed TO A 
albedos by constraining the surface albedo and in this case, it did lead systematically to an increase of 
32 solar irradiance values in the continental US. Other potential sources of the bias were examined. 

Differences in elevations of station locations and of mean elevations of NASA/POWER grid cells 
34 appeared at best to explain only a small portion of the bias (Fig. 5B), and elevation differences also 
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appeared to have little influence on correlations between SRAD 0 b and SRAD NP (Fig. 5A). Another 
2 possible explanation for the bias concerns locations of the weather stations. Stations associated with 

airports or agricultural research centers may have been located in open areas with a clear field of view and 
4 low probability of cloud cover, while NASA/POWER grid cells may have included mountainous areas 
that experienced lower daily radiation due to greater cloud cover. Reflection from clouds can increase 
6 irradiance on a scale of minutes, but such effects tend to be cancelled out when the sun is obscured by 
clouds (Pfister et al., 2003) and thus seem unlikely to contribute to the bias. The occurrence of sporadic 
8 excessively high values of SRAD 0 b also might have biased the mean value of SRAD 0 b, or it may 
evidence a tendency of some automated instruments to overestimate SRAD. As a partial test for this 
10 problem, the mean of SRAD 0 b was re-calculated after limiting all values of SRAD to a maximum of 0.8 
of Qo. The resulting mean was 17.25 MJm 2 d _I as compared to 17.32 MJm" 2 d 1 for the dataset as used in 
12 the rest of the paper, where values greater than Q 0 were simply excluded. Thus, we suggest that errors in 
the prescription of the background aerosol for the GEWEX SRB solar irradiance v2.81 are the largest 
14 contributor to the noted systematic bias. 

The NASA/POWER data accurately reproduced the variability in solar radiation data, did not 
16 show major variation related to effects of elevation, and are readily available via the Internet for a time 
span from 1984 onward. They thus show excellent potential as a source of solar radiation data for diverse 
18 applications. However, the differences in mean values of SRAD 0 b and SRAD NP were a concern, but 

appear to be at least partially addressed in the GEWEX SRB solar v3.0 that is now available. If analysis 
20 of the new version largely explains the differences noted in this paper, the NASA/POWER data would 
appear to be superior to data from weather generators. 

22 A previous paper comparing NASA/POWER daily temperature data to COOP data found larger 

differences between these two sources (White et al., 2008). Several hypotheses can be forwarded to 
24 explain why the NASA/POWER solar data may show greater consistency than the temperature data. The 
first and foremost is that the estimate of solar radiation is far more dependent upon the structure and 
26 variability of the actual cloud fields than on accuracy in the underlying meteorological fields from the 

atmospheric assimilation data sets. Those cloud fields are directly observed by the satellite measurements. 
28 In fact, only the total water vapor profile is required for the solar radiation calculation and even 10-30% 
errors in the water vapor profile amount (which are probably correlated with surface temperature errors in 
30 the assimilation) correspond to solar radiation errors < 1%. Time of observation bias should have been 
less of a problem since observed solar radiation data presumably were based on a midnight to midnight 
32 integration period. Instrument siting errors may have less effect on solar radiation data than on 

temperature data. We emphasize that the algorithms used to create the datasets are subject to periodic 
34 review and improvement as noted by the effect of the improved aerosol inputs described above. 
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Additional improvements in the satellite calibration, sampling size and solar algorithm is anticipated and 
2 in 2012, it is planned to reduce the grid cell size to 0.5°. 

In working with datasets from automatic weather stations, accessed via the Internet, numerous 
4 problems with data quality were encountered. Together, these reinforce concerns over the management of 
daily weather data (e.g., Davy and Pielke, 2005; Holder et al, 2006; Pielke et al., 2007). Various data 
6 checking procedures exist, but they do not appear to be used by all data providers. The variable quality 
control provides another argument in favor of using a single, well-documented data source such as offered 
8 by NASA/POWER. 

10 5. Conclusions 

12 Considering the constraints inherent with its coarse grid size of 1° x 1° of latitude and longitude, 

the NASA/POWER solar radiation data compare favorably with data reported from automatic weather 
14 stations. In terms of representing historic variability on time scales of a few days, they appeared superior 
to values estimated using the WGENR weather generator. However, the means of the SRAD NP data were 
16 often 1 to 2 MJm 2 d"’ lower than SRAD 0 b, and this discrepancy merits further investigation. 

NASA/POWER data are available for over 25 years with global coverage and are continuously 
1 8 being updated and improved. They represent a valuable source of solar radiation data for research 
concerned with regional to global geographic scales. 

20 
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Captions for Figures 

2 Figure 1. Comparisons of correlations and differences between mean values of daily solar radiation data 
from different sources. A. Correlation between SRAD 0 b and SRAD NP . B. Correlation between 
4 SRADob and SRADwg- C. Difference between SRAD 0 b and SRAD NP . D. Difference between 

SRADob and SRADwg- 

6 

Figure 2. Comparisons of solar radiation values (in units of MJm 2 d _1 ) from the various sources of SRAD 
8 data for Immokalee, Florida from 1 Jan. 1998 to 31 Dec. 2004. The diagonal lines represent a 1:1 

relation. A. SRADob vs SRADn P . B. SRADob vs SRADwg- C. SRADob vs SRADco- D. SRADco 
10 vs SRAD np . All values of r 2 are significant at the P < 0.001 level. 

12 Figure 3. Root mean square errors (RMSE) for prediction of observed solar radiation from the other 

sources. A. SRADob vs SRADn P . B. SRADob vs SRADwg- C. SRADob vs SRADco- D. SRADob 
14 vs Qo. 

16 Figure 4. Density plots comparing values of solar radiation (MJ m 2 d~') for 295 stations in the continental 
US. A. Weather station vs. NASA/POWER. B. Weather station vs. WGENR-generated. Count 
18 ranges are 1 = 1 to 10 paired values; 2 = 1 1 to 100; 3 = 101 to 1000; 4 = 1001 to 10,000. 

20 Figure 5. Annual variation in solar radiation based on seven-day averages. A. Daily means across all 
locations for SRADob, SRAD np and SRADwg- B. Mean difference between SRADob and 
22 SRAD np for five latitude bands from less than 30°N to greater than 45°N. C. Relative error 

[(SRAD np - SRAD 0 b)/SRADob] for the five latitude bands. 

24 

Figure 6. Relation between difference in elevation of weather station and of the corresponding grid cell 
26 for the NASA/POWER dataset and two indicators of reliability for all stations. A. Correlation 

between SRADob and SRAD NP . B. Difference between mean of SRADob and mean of SRAD np . 

28 

Figure 7. Comparison of solar radiation data from NASA/POWER and estimated with WGEN using daily 
30 weather data from NOAA COOP stations. A. Correlation between the two sources for each grid 

cell. All correlations are significant at the P < .001 level. B. Mean difference between the NOAA 
32 COOP and NASA/POWER sources. 
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Supplement Fig. 1. Variation in solar radiation for five stations where the data showed patterns suggesting 
2 problems in instrumentation. A. Semmes, AL. B. Avondale, CO. C. Novelty, MO. D. 

Calipatria/Mulberry, CA. E. Manteca, CA. 

4 

Supplement Fig. 2. Variation in the ratio of observed solar radiation to extraterrestrial (Q 0 ) for five 
6 stations where the data showed patterns suggesting problems in instrumentation. The horizontal 

reference line at 0.8 indicates an approximate upper limit for any location. A. Semmes, AL. B. 

8 Avondale, CO. C. Novelty, MO. D. Calipatria/Mulberry, CA. E. Manteca, CA. 

10 
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Table 1. 


2 Daily meteorological variables available on a global 1° grid through the NASA/POWER project. 


Variable 

Source 

Time span 

Availability 
from present 
date 

Daily 

Maximum and 
minimum 
temperatures, 
Daily averaged 
temperature 

Goddard Earth Observing System 
(GEOS) assimilation model version 4 

Goddard Earth Observing System 
(GEOS) assimilation model version 5 

January, 1983 to December 
2007 

January 2008 to present 

online 
< 1 week b 

Precipitation 

Satellite & ground observations from 
the Global Precipitation Climatology P 
project (GPCP) 

January, 1997 to present 1 * 

< 2 months 

Solar radiation 

Satellite observations 
GEWEX SRB v3.0* 

FLASHFlux 

July 1983 to 
December 2007 
January 2008 to Present 

online 
< 1 week b 

Dewpoint 

temperature 

Goddard Earth Observing System 
(GEOS) assimilation model version 4 

January, 1983 to December 
2007 

online 


Goddard Earth Observing System 
(GEOS) assimilation model version 5 

January 2008 to present 

< 1 week b 


4 a The most current GPCP file is August 2009. 
b Data files are updated daily. 

6 * v3.0 came available after the analysis presented here was completed; v2.81 corresponding to White et 

al., (2008) was used for this study. 
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Table 2. 

2 Mean, minimum and maximum values of solar radiation for the four solar data sources. The 

NASA/POWER and COOP data are for 855 locations on a 1° latitude and longitude grid covering 
4 the continental US and representing a time series from 1983 through 2005. Elevations correspond 

to mean values of grid cells for the NASA/POWER dataset and to reported values for COOP 
6 stations. 


Data source 

Mean 

Minimum 

Maximum 

Sites with automated stations (N = 295) 

Automated stations 

1 7.4 a 

0.2 a 

43. 0 a 

NASA/POWER 

16.2 

0.3 

34.1 

Generated based on station data 

16.5 

1.3 

33.2 

Generated based on COOP data 

16.8 

1.3 

33.2 

Difference between NASA/POWER and Automated 
Stations 

-1.2 

-30.4 

29.6 

Elevation difference (m): NASA/POWER -Automated 
Stations 

+219 

-847 

368 

Entire US based on COOP stations (N = 855) 

NASA/POWER 

15.0 

0.1 

34.2 

Generated based on COOP data 

15.9 

1.3 

33.2 

Difference between NASA/POWER and COOP-based data 

-0.9 

-31.1 

26.1 

Elevation difference (m): NASA/POWER - COOP 

85 

-1580 

1270 


8 a Minimum and maximum values of SRAD from automated stations reflect minimum and maximum 
values imposed in processing data. 
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